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The Derivation of the Exact Internal Energies for Spin Glass 
Models by Applying the Gauge Theory to the Fortuin-Kasteleyn 

Representation 

Chiaki Yamaguchi 
Kosugichou 1-359, Kawasaki 211-0063, Japan 

We derive the exact internal energies and the rigorous upper bounds of specific heats for 
several spin glass models by applying the gauge theory to the Fortuin-Kasteleyn representa- 
tion which is a representation based on a percolation picture for spin-spin correlation. The 
results are derived on the Nishimori lines which are special lines on the phase diagrams. As 
the spin glass models, the ±J Ising model and a Potts gauge glass model are studied. The 
present solutions agree with the previous solutions. The derivation of the solutions by the 
present method must be useful for understanding the relationship between the percolation 
picture for spin-spin correlation and the physical quantities on the Nishimori line. 

§1. Introduction 

The theoretical studies of spin glasses have been widely doneP There are special 

lines on the phase diagrams for several spin glass models, where the lines are called 

the Nishimori line.^®® It is known that several physical quantities and several 

bounds for physical quantities are exactly calculated on the Nishimori line by using 

gauge transformations. The exact internal energies and the rigorous upper bounds of 

specific heats for several spin glass models have already been derived on the Nishimori 
HneJlM The 

aim of this article is to derive the exact internal energies and the 
rigorous upper bounds of specific heats for several spin glass models by applying the 
gauge theory to the Fortuin-Kasteleyn representation. As the spin glass models, the 
±J Ising modeP 1 '^ and a Potts gauge glass modeP 1 are studied. The Potts gauge 
glass model is a more complex version of the ± J Ising model. 

The Fortuin-Kasteleyn representation is a representation based on a percolation 
picture for spin-spin correlation.^ 'EJ 1 ^ By using the Fortuin-Kasteleyn represen- 
tation, spin-spin correlations are directly treated. In the previous methods, the 
solutions have been directly calculated from the Boltzmann factorP ™™ Instead, in 
the present method, the solutions are directly calculated from the Fortuin-Kasteleyn 
representation. 

In this article, gauge transformations are used. The gauge transformations are 
treated in Refs. .15]) .I5|) .[TU|) for example. It is known that a gauge transfor- 

mation has no effect on thermodynamic quantities P 1 

This article is organized as follows. The Fortuin-Kasteleyn representation is 
briefly explained in £j2j The solutions for the ±J Ising model are obtained in $31 
The solutions for a Potts gauge glass model are obtained in JJH This article is 
summarized in £J5J Appendices A and B are attached in order to make this article 
self-contained. 
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§2. The Fortuin-Kasteleyn representation 

We briefly explain the Fortuin-Kasteleyn representationF^® The Fortuin-Kasteleyn 
representation introduces auxiliary variables called graph G. The graph G is a state 
by the weights between spins which are directly connected by the interaction. The 
partition function Z is expressed in the double summation over state S and graph 
GaP 

Z = J2V(G)A(S,G), (2-1) 

S,G 

where A(S,G) is a function that takes the value one when S is compatible to G 
and takes the value zero otherwise. A bond that probabilistically connects two spins 
by the weight of graph is especially called the active bond. A graph consists of a 
set of active bonds. The active bond is fictitious, and is used in order to generate 
a cluster composed of spins, which is often referred to as the Fortuin-Kasteleyn 
cluster. V(G) is the weight for the graph G. The partition function Z is expressed 
as J2g W(G), where W(G) = V(G) Y^s ^K^' Gf). This partition function is expressed 
in the summation over graph G instead of state S. Since the weight for the graph G 
is used, the study of the percolation problem of the graph G is expected to have a 
physical significance. This representation for graph G is called the Fortuin-Kasteleyn 
representation. 

We define the number of active bonds as n&. The number of states for the active 
bond number, i?(n&), is given by 

{G\n b (G)=n b } S 

By using the f2(rib), the partition function is expressed as 

N B 

Z=Y, n(n b )V(n b ) , (2-3) 

"6=0 

where V(n b ) is the weight for the active bond number, and Nb is the number of 
nearest-neighbor pairs in the whole system. 

If two spins are on the same cluster, the two spins are correlated. If two spins 
are not on the same cluster, the two spins are not correlated. In the ferromagnetic 
Ising model, the percolation transition point of the Fortuin-Kasteleyn cluster agrees 
with the phase transition point On the other hand, in the ±J Ising model, the 
percolation transition point of the Fortuin-Kasteleyn cluster disagrees with the phase 
transition pointP^ 1 Instead, it is pointed out that, in the ±J Ising model, there is 
a possibility that the percolation transition point of the Fortuin-Kasteleyn cluster 
agrees with a dynamical transition point For the applications of the Fortuin- 
Kasteleyn representation, the Swendsen-Wang algorithm^ jg probably the prime 
example. This algorithm is a Markov chain Monte Carlo method. By performing 
this algorithm, the Fortuin-Kasteleyn clusters are generated, and the states on each 
cluster are simultaneously updated. This algorithm produces a faster thermal equi- 
libration when this algorithm is applied to the ferromagnetic Ising modelP^ In this 
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article, we concentrate ourselves on the number of the active bonds, which generate 
the Fortuin-Kasteleyn clusters, and the fluctuation of the number of active bonds. 

§3. The ±J Ising model and the present results 

The Hamiltonian for the ±J Ising model, H, is given bjP®^ 

H = — Jij^iSj ' (3d) 
(hi) 

where denotes nearest-neighbor pairs, Si is a state of the spin at the site i, 

and Si = ±1. Jij is a strength of the exchange interaction between the spins at the 
sites i and j. The value of Jjj is given with a distribution P(Jjj). The distribution 
P(Jij) is given by 

P(J id ) = P 6j iJ ,j + (1 - p) 6j iJt -j , (3-2) 

where J > 0, and 5 is the Kronecker delta, p is the probability that the interaction is 
ferromagnetic, and 1 — p is the probability that the interaction is antiferromagnetic. 
By using Eq. flHJgp , the distribution P(Jij) is written asP®'® 

p(J - ) = 2«7)' J - = ±J ' < M > 

where /3p is given by®'®® 

Pp-ij^. (3-4) 

When the value of /3p is consistent with the value of the inverse temperature f3, the 
line on the phase diagram for the temperature T and p, where Eq. (|3-4p is satisfied, 
is called the Nishimori line. 

A gauge transformation® 'IS 1 'El 1 given by 

OiOj , Si ->■ SiOi (3-5) 

is used where o~i = ±1. By using the gauge transformation, the Hamiltonian % part 
becomes % — > %, and the distribution P(Jjj) part becomes 



II p ( j mO 



[2cosh(/W)]^ 

2^[2cosh(/3 P J)]^ ' 1 J 

where iV is the number of sites. 

For the ± J Ising model, V{n^) is given by 

V(n b ) = (e 2/3J - \) n b e -N B l3J ; (3 . 7) 

where j3 = T is the temperature, and ks is the Boltzmann constant. The 

way of deriving Eq. (|3-7|) is described in Appendix A. We define the probability 
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for putting the active bond as -Pfk- The value of Pfk depends on the exchange 
interaction and the states of spins.®®" 1 ^'®'^ For the ±J Ising model, Pfk is 
given bfEUSl 

P FK (S U Sj, Jij) = 1 - e-fUuSiSj-PlJul (3 . 8) 

The way of deriving Eq. (j3-8|) is also described in Appendix A. By using the gauge 
transformation, the Pfk part becomes Pfk — > Pfk- 
The internal energy E is given by 

E = ~\hxZ\ R , (3-9) 

where [] R denotes the random configuration average. By using Eqs. (j2-3|) . ()3-7p and 
(|3-9p . we obtain 

2 Je 2 P J 

E = N B J- eW _ i [(n b ) T ] R , (3-10) 
where ()r denotes the thermal average. [(n^rlp is given by 

[(n b ) T ] R = J2l( p FK(Si, Sj, JijMn - (3-11) 

When (3 = /3p, [(Pfk(£ij Ji,j))T\R is obtained by using the gauge transformation 



{J ; , m } (i, m > C 
1 



J) 



" 2^[2coBh(/3 P J)]^ 

= tanh(/W) , (3-12) 

where /3p = l/kpTp, Tp is the temperature on the Nishimori line. By using 
Eqs. ([HTTP and (|3 : T2l) . we obtain 

[(n b ) T } R = N B t<mh(f3 P J). (3-13) 

By using Eqs. (|3T(jp and (|3T3|) . the internal energy E is obtained as 

E = - N B J 'tanh(/9p J) . (3-14) 

This solution is exact, and is equivalent to the solution in Ref. [2]). 
The specific heat C is given by 



( '- hjp 2 -^[lvZ) l; . CM5) 
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By using Eqs. ([23]), (g£fl) and ([345]), we obtain 

C = k B (f3J) 2 cosech 2 (/3J){-[(n b ) T ] R 

+ e 2PJ ([(nl) T } R -[(n b ) 2 T ] R )}. (3-16) 

K"|Hi? is given by 
[{nl) T ] R 

= ^ ^jKffrKpSj) Sj, Jjj)PFK(Sk, Si, Jfc,f)(l — &i,k$j,l — 5i,l0~j,k) 
(hi) {k,l) 

+ Pfk(Si, Sj, Jij){5i,kSj,i + $i,i$j,k))T]R ■ (3-17) 

By performing a similar calculation with the calculation in Eq. (|3-12p . we obtain 

l(n 2 } T ] R = N B (N B - l)tanh 2 (/3 P J) 

+ iV B tanh(/3 P J). (3-18) 

By applying the Cauchy-Schwarz inequality, we obtain 

l(n b ) 2 ] R > [(n b )] 2 R = N 2 B tanh 2 (/3 P J) . (3-19) 

Therefore, by using Eqs. ()3T3p . ()3T6p . (|3T8p and (|3-19p . we obtain the upper bound 
of the specific heat C as 

C < k B N B (/3pJ) 2 sech 2 (PpJ) . (3-20) 
This solution is rigorous, and is equivalent to the solution in Refs. E]),[3]). 

§4. A Potts gauge glass model and the present results 

The Hamiltonian for a Potts gauge glass model, T~L, is given by® 



H = " Z) E e^ ( "« + *-* )rw , (4-1) 



where % is a state of the spin at the site i, and qi = 0, 1, . . . , q — 1. z/j j is a variable 
related to the strength of the exchange interaction between the spins at the sites i 
and j, and f = 0, 1, . . . , q— 1. g is the total number of states that a spin takes. The 
value of z/jj is given with a distribution P(fjj). The distribution P(Vjj) is given by 

= P^,o + (1 - <^,o) • (4-2) 

The normalization of P{vij) is given by 

]Tp(^) = l. (4-3) 
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When Ui t j = for all pairs, the model becomes the ferromagnetic Potts model. 
When q = 2, the model becomes the ±J Ising model. By using Eqs. (|4-2p and (|4-3p . 
the distribution P(i>ij) is written as®'^ 



P(u itj ) = Ae 
where A and (3p are given bj^'^ 



q Z^r 4 =1 \"l>,]> 



A 



e^ (9_1) + (q- l)e 



1 



1 



(4-4) 



(4-5) 



(4-6) 



respectively. When the value of /3p is consistent with the value of the inverse temper- 
ature /?, the line on the phase diagram for the temperature T and p, where Eq. (|4-6p 
is satisfied, is called the Nishimori line. 

We use representations: Aj = e « and = Je i M 1,3 . A gauge transfor- 

mation®'^ given by 



j( r i,j) v j( r i,j).. C l r i,3 ,, r i,3 

J i,j ~? J i.i H 



1,3 



(4-7) 



is used where //j = e" 9 ', is an arbitrary value for the spin state at the site i, 
and q~i = 0, 1, . . . , q — 1. By using the gauge transformation, the Hamiltonian Ti part 
becomes % — > and the distribution P(vij) part becomes 



A NB e~ ^ r ».i =1 ^ 
A n b 



i,3 1 



dp 



For the Potts gauge glass model, T^(nb) is given by 

y(n fe ) = (e^ - l) n *e ~ . 



(4-8) 



(4-9) 



The way of deriving Eq. (|4-9p is described in Appendix B. For the Potts gauge glass 
model, Pfk is given bjfl^ 



Ti,j=l 



(4-10) 



The way of deriving Eq. (|4-10p is also described in Appendix B. By using the gauge 
transformation, the Pfk part becomes Pfk - ► Pfk- 

By using Eqs. (|2-3|) . (|3-9|) and (|4-9j) . the internal energy E is given by 



(4-11) 
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[Wt]h is given by 

,Qj,Vi,j))T]ll- (4-12) 

When = f3 P , [(P FK (qi ' 1j' u i,j))T\R is obtained by using the gauge transformation 
aP* 

{(PFK(qi,qj,Vi,j))T]R 



2^ 11 ^W,mJ e -/3p«({</«},K m }) 
Km} (i,m) Mft} c 

/iATh _ 



Km} {?i} 

= -^1 7> ( 4 ' 13 ) 

where f3p is the inverse temperature on the Nishimori line. By using Eqs. (|4-12[) and 
(|4- 13j) . we obtain 

By using Eqs. (|4-lip and ()4-14p . the internal energy E is obtained as 

jVgJ_ NsJeM 
q ePp J + q-l y J 

This solution is exact, and is equivalent to the solution in Ref. H|). 

By using Eqs. (|2-3|) . (|3- 15|) and (|4-9|) . the specific heat C is given by 

C = kB {^f) cosech2 f^~V _ ^ n6 )T]iJ 

+ ^ J ([{nl) T ) R -[{n b )l) R )}. (4-16) 

[( n b)T]i? is given by 

= ^2^2[{PFK{qi,qj,Vi,j)PFK{qk,qi,Vk,i){l ~ $i,k°~j,l ~ h,io~j,k) 
(id) (k,l) 

+ PFK{qi,qj,Vi,j){o~i,kO~j,i + S it i5 jjk )} T ] R . (4-17) 
By performing a similar calculation with the calculation in Eq. ()4-13[) . we obtain 

N B (N B - l)(e^ J - l) 2 



[{<) 



T R 



{e?r J + q-l) 2 

+ yf J -'>. (« 8) 
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By applying the Cauchy-Schwarz inequality, we obtain 

N 2 ( J3 P J _ i\2 

l(n b ) 2 ] R > [(n 6 )]| = { f pJ - - _ l )2 ■ (4-19) 

Therefore, by using Eqs. (|4-14p . (|4-16p . (|4- 18j) and (|4-19p . we obtain the upper bound 
of the specific heat C as 

(ePP J + q - l) 2 

This solution is rigorous, and is equivalent to the solution in Ref. H|). 



§5. Summary 



We derived the exact internal energies and the rigorous upper bounds of spe- 
cific heats for the ±J Ising model and a Potts gauge glass model by applying the 
gauge theory to the Fortuin-Kasteleyn representation. The results were derived on 
the Nishimori lines. The present solutions agreed with the previous solutions in 
RefsP'^'Sl 1 The Fortuin-Kasteleyn representation is a representation based on a 
percolation picture for spin-spin correlation. The derivation of the solutions by the 
present method must be useful for understanding the relationship between the perco- 
lation picture for spin-spin correlation and the physical quantities on the Nishimori 
line. 



Appendix A 

The weight and the probability for active bond in the ± J Ising model 

We will derive Eqs. (|3-7p and (|3-8p . The framework for the way to derive 
Eqs. (|3-7|) and (|3-8|) is described in Ref. [7J). We define the weight of two spins 
as w(Si,Sj,Jij). w(Si,Sj,Ji tj ) is given by 

w(Si, Sj, J it j) = exTp({3J itj SiSj) . (A-l) 

We define the weight for JijSiSj = J as u> para - We obtain 

w P ara(5'i, Sj, J{,j) = exp(/3 J) . (A-2) 

We define the weight for J^jSiSj = — J as w a nti- We obtain 

uwOSi, Sj, Ji,j) = exp(-/3J) . (A-3) 

We define the weight of graph for connecting two spins as w(g conn ). We define the 
weight of graph for disconnecting two spins as u'(gdisc)- We are able to write 

w para (S'i, Sj, Jij) = w(g conn ) + w(g d isc) , (A-4) 
Wanti(Si, Sj, J^) = w(g disc ) . (A-5) 
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By using Eqs. (J£j2j), flEI , <lA4l) and (jA^j) . we obtain 

w(g conn ) = exp(/3J) - exp(-/3J) , (A-6) 
w{g d i sc ) = exp(-/3J) . (A- 7) 

By using Eqs. (|A-6p and (|A-7p . we obtain the weight V(nb) for the active bond 
number n& as 

y(n b ) = (eP J - e~ pj ) nb { e ~P J ) N B- n i- (A-8) 

The above equation is equal to Eq. ()3-7p . We define the probability of connecting 
two spins for JijSiSj = J as -P pa ra (^conn)- We define the probability of connecting 
two spins for JijSiSj = — J as -Panti (strain)- We are able to write 



w{g c , 

w(g 

cram ) + ^(Sdisc) ' 
-Panti(Sconn) = 0. (A- 10) 



-'para vSconn ) — „../"„ \ , „..t „ \ ' V^'^J 



By using Eqs. (TM . (jA7ll . (iA9l) . ([AdOi we derive Eq. (^81) . 

Appendix B 

TTie weight and the probability for active bond in a Potts gauge glass model 



We will derive Eqs. (|4-9p and (|4-lUp . The framework for the way to derive 
Eqs. ()4-9p and (|4-10p is described in Ref. [7J). We define the weight of two spins as 
w(qi,qj,Vij). w{qi,qj,Uij) is given by 



w(qi,qj,Vi d 

9-1 



exp 



cxp 



2m 

— I "t j + Qi-Qj 



We define the weight for u^j + qi — qj = as w pa ra- We obtain 

w pai r a {qi,qj,v i:j ) = exp 

We define the weight for v^j + q% — qj ^ as w a nti- We obtain 

w ant i{qi, qj,Vi,j) = exp I — — 



(B-l) 



(B-2) 



(B-3) 



We define the weight of graph for connecting two spins as w(g COTin ). We define the 
weight of graph for disconnecting two spins as w(sdisc)- We are able to write 



^para(%i qj > — ^(Sconn ) + ^(Sdisc) , 

Wantife, Qj, = ^(ffdisc) ■ 



(B-4) 

(B-5) 
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By using Eqs. (|B^2|> . fETS]) . (jB~4)l and (jB~5l> . we obtain 



w(g conn ) = exp 



9 



™(#disc) = eXP (^~^J ■ ( B ' 7 ) 

By using Eqs. (|B-6p and ()B-7p . we obtain the weight V(rib) for the active bond 
number n\, as 

PJ(q-l) 0J 0J „ 

y(n 6 ) = [e^^~ - e"p ( e — r)^-"* _ ( B . g ) 



The above equation is equal to Eq. (|4-9p . We define the probability of connecting 
two spins for Vij + qi — qj = as -P P ara(<?conn)- We define the probability of connecting 
two spins for i/y + qi — qj ^ as -P a nti(5conn)- We are able to write 



w(g C onn) 

w(g 

conn ) + ^(^disc) ' 
Pantifeonn) =0. (B-10) 



-'para \yconn J — „../"„ \ i „../"„ \ ' \D-if) 



By using Eqs. (uTBD . (jB~7l) . (jB~9]l and p40]) . we derive Eq. (jCTD . 
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